Enabling technology for global 3D <b>+</b> 3V hybrid-Vlasov simulations of near-Earth space

We present methods and algorithms that allow the Vlasiator code to run global, three-dimensional hybrid-Vlasov simulations of Earth's entire magnetosphere. The key ingredients make Vlasov at magnetospheric scales possible are sparse velocity space implementation spatial adaptive mesh refinement. outline algorithmic improvement semi-Lagrangian solver for six-dimensional phase quantities, discuss coupling Maxwell equations' solvers in a refined mesh, provide performance figures from simulation test runs demonstrate scalability this system full runs.